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INTRODUCTION 


The recent discovery of asymptotically safe quantum field theories in four space-time 
dimensions offers new directions for both fundamental theory and model building [T]. The 
main novelty of these theories is the presence of an exact interacting ultraviolet (UV) fixed 
point. This distinguishes them from asymptotically free theories where the UV fixed point 
is non-interacting Hi- Also, no additional symmetry principles such as space-time su¬ 
persymmetry |1] are required to ensure well-defined and predictive theories in the UV [5]. 
Instead, the fixed point arises dynamically through renormalisable interactions between non- 
Abelian gauge fields, fermions, and scalars, and in a regime where asymptotic freedom is 
absent. Crucially, the potentially dangerous growth of the gauge coupling towards the UV 
is countered by Yukawa interactions, while the notorious growth of Yukawa and scalar cou¬ 
plings is tamed by the fluctuations of gauge and fermion fields. This has led to theories 
with “complete asymptotic safety”, meaning interacting UV fixed points in all couplings [1] . 
This is quite distinct from the more conventional setup of “complete asymptotic freedom” 
M, where the UV dynamics of Yukawa and scalar interactions is brought under control 
by asymptotically free gauge helds; see pen] for recent studies. 

In order to qualify as fundamental theories, UV hxed points must reside in the “physical” 
regime of parameter space where the vacuum is stable quantum-mechanically. The central 
question we wish to address in this paper is: Does the quantum vacuum remain stable at an 
interacting UV hxed point? It is well-known that huctuations can induce the spontaneous 
breaking of symmetry d la Coleman and Weinberg m- We will answer the question to 
the positive by performing a detailed analysis of the theory’s phase diagram. We also offer 
results for key features of the theory including characteristic energy scales, the running 
couplings, the separatrix, and the resummed Coleman-Weinberg potential. 

The rest of the paper is organised as follows. We recall the main features of asymptotically 
safe gauge-Yukawa theories, which is followed by a detailed study of the phase diagram and 
the UV completion of the theory including analytical expressions for running couplings 
along UV-safe trajectories (Sec. 0. We then perform a stability analysis of the scalar field 
potential, both classically and quantum-mechanically, also offering analytical and numerical 


results for the resummed Coleman-Weinberg effective potential (Sec. Ill), followed by our 


conclusions (Sec. IV). 


II. PHASE DIAGRAM AND UV COMPLETION 

In this section, after reviewing the set-up and salient perturbative properties of asymptot¬ 
ically safe gauge-Yukawa theories introduced in [I] , we discuss the phase diagram including 
the RG evolution of the quartic self-couplings. We also derive explicit expressions for the 
UV-safe trajectories away from the UV stable fixed point including along the UV-IR con¬ 
necting separatrices, and further quantities required for the subsequent stability analysis. 
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A. Preliminaries 


Following pp, we consider a massless quantum field theory with SU{Nc) gauge helds 
and held strength (a = 1, • • • , N^, — 1), Np havors of fermions Qt {i = 1, ■■■ , Np) in 
the fundamental representation, and a. Np x Np complex matrix scalar held H invariant 
under U{Np)p x U{Np)p rotations and uncharged under the gauge group. The fundamental 
Lagrangian is given by the sum of the Yang-Mills term, the fermion and scalar kinetic terms, 
the Yukawa interaction, and scalar self-interaction terms. 


L = - f Tr + Tr (Q tip Q) + Tr {d^H^ d^H) 

+y Tr {QpHQp + QrH^Ql) - uTt - v (Tr H^Hf , 


( 2 . 1 ) 


where the decomposition Q = Qp -|- Qp with Qp/p = |(1 ± 75 )Q is understood. The trace 
Tr indicates the trace over both color and havor indices. The model has four classically 
marginal coupling constants in four space-time dimensions given by the gauge coupling g, 
the Yukawa coupling y, the quartic scalar couplings h and the ‘double-trace’ scalar coupling 
V, which we write as 


ag — 


g^Nc 

(47r)2 


ay — 


y^Nc 

(47r)2 


Oih = 


u Np 
(47r)- 


CXt) — 


vN^p 

(47r)^ 


( 2 . 2 ) 


We have also added the appropriate powers of Nc and Np in the normalization of the 
couplings. The shorthand notation (3i = dtai with i = (^f, y, h, v) is employed to indicate the 
/^-functions for the couplings (2.2). In the Veneziano limit where both Nc and Np are large 
but their ratio hxed, the parameter 


Np 


e = 


N 


c 


11 

y 


(2.3) 


can take any real value. In the large-Y limit, the perturbative renormalisation group equa¬ 
tions for the couplings (2.2) have been obtained in [12] in dimensional regularisation, also 


using the results HMH). In terms of (|2.3[) they are given by 


26 


y ~ + A 25 -|- —e — 2 — e 


11 




a. 


701 53 


112 


27 






——e ) Q!g —^(11 -h 2e) agay + -(11 -h 2e) (20 3e)a 

-|- 2e) ay — 6 a^ 

20e - 93 
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■i 2 7. . ^ / 385 23 

Oq + (49 + 8e)ag% - I — + — e + — 


(2.4) 


a. 


(2.5) 


ay — (44 -|- 86)ayafi 4q;^ 


(5h — “(11 + 2e) tty -|- Aah{ay 2ah ), 
jS-u = 12cr^ -|- 4q:^ [a^ -)- 4cr/i -|- ay) . 


( 2 . 6 ) 

(2.7) 


for (3g,(3y,(5h and jiy up to (3,2,1, l)-Ioop order, respectively. In the terminology of P we 
refer to this as the next-to-next-to-Ieading order (NNLO) approximation. The NLO approx¬ 
imation corresponds to the approximation in which the (2,1, 0, 0)-loop terms for f3g, (3y, (3h 
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coupling 

order 

in perturbation theory 

ag 

1 

2 

2 

3 

Oy 

0 

1 

1 

2 

Oih 

0 

0 

1 

1 

Oiy 

0 

0 

1 

1 

approximation level 

LO 

NLO 

NLO' 

NNLO 


Table I: Relation between approximation level and the loop order to which couplings are retained 
in perturbation theory. 


and are retained. As discnssed in [T] , this ordering of pertnrbation theory is also favonred 
by Weyl consistency conditions [laiisiiiz]. Additionally, we consider the NLO' approxima¬ 
tion in which the (2,1,1, l)-loop terms are retained (see Tab. [T]). We stress that the dynamics 
of the scalar conplings, central for our study of vacuum stability, arises for the hrst time 
within the NLO' and NNLO approximations. 


The theory (2.1) is renormalisable within perturbation theory. For e < 0 it also displays 


asymptotic freedom in the gauge sector. When e is positive, asymptotic freedom is lost and 


the gauge sector becomes QED-like. In [T], it has been established that the theory ( 2 . 1 ) 
develops an exact, interacting UV hxed point in all four couplings, strictly controlled by 


perturbation theory, provided the parameter (2.3) is chosen to be 


0 < e < 1 . 


(2.8) 


The existence of an interacting UV hxed point ensures that the theory (2.1) remains predic 


tive to highest energies at the hxed point. Another feature of this theory is that the scalar 
sector avoids a triviality bound due to residual interactions in the UV, meaning that the 
scalars can be viewed as elementary. 


B. Fixed points 


In |1], it has been shown that the gauge-matter system (2.4) - (2.7) may display three 


diherent ultraviolet hxed points in the regime (2.8). In an expansion in the small parameter 
e, the hrst hxed point has the coordinates 


_ 26 I 23(75245-130687^) ^2 
“9 “ 57 ^ 


370386 


e- + (!l(e3) 


Of, = 


19 

V ^-1 


19 


f 43549 
I 20577 


e + 0(e^), 


23007 ^ 

6859 


) + 0((^) 


(2.9) 


^vl 19 


^( 2 ^/^ - \/20 + 6 ^) e + C>(e 2 ). 


We refer to the hxed point (2.9) as FPi. There are two more hxed points in the perturbative 


domain (2.8), a second fully interacting hxed point FP 2 and a partial hxed point FP 3 . The 























5 


coordinates of FP 2 differ from FPi in (2.9) in the coordinate whose hxed point becomes 


^v2 19 


(IspS, + V 20 + 6y23)e + O(e^) 


( 2 . 10 ) 


instead. The second hxed point comes about because the double trace coupling does not 

t-u to all 


couple to the RG how of {ag, ay, ah), and its own RG how (2.7) is quadratic in a^ 
orders in l/Np and 1/Nc- Therefore, a^ may even interpolate between these two hxed 
points without ahecting the hxed point for (a^, ay, ah)- 

Finally, a partial hxed point FP 3 arises where all but the double-trace coupling ay settle 
onto a hxed point within the perturbative domain. The double-trace coupling continues to 
evolve logarithmically towards a perturbative Landau pole. Furthermore, ah settles at a 


negative hxed point value while a* and a* remain equal to their values at (2.9) to this order 


m e. 


Towards UV momentum scales, the partial hxed point FP 3 is then characterised by 


a 


h2 


= -{V^+l) ^ + 0(e 


19 


Civin) = ay(fio) + - K tan ( neln — 

4 V ho 


( 2 . 11 ) 


where k = ^^/6y/2K--^ ^ 0.624. We have chosen the scale po as a free parameter with 


otv{l^o) = (2ev^)/19 ~ 0.505 e. The running for ay in (2.11) is an exact solution of its 


beta function (2.7) provided that a„, and ah have settled on FP 3 to leading order in e. 


From (2.11) we observe that ay reaches its perturbative UV Landau pole at 


vr 


ALandau ~ hO exp -- > /iQ . 

2Ke 


( 2 . 12 ) 


A similar Landau pole is reached towards low energies. 


C. Phase diagram 

In [1], the phase diagram of the model has been discussed at NLO accuracy. Here we 
extend the study to NLO' and NNLO order including ehects due to the running scalar 
couplings, and illustrate results in Figs. [T}[^ and which show portions of the phase diagram 
in the vicinity of the UV and IR hxed points. 

Specihcally, Fig. shows the RG trajectories for the couplings {ag, ay, ah), including 
the UV and IR hxed points and the UV-complete trajectory (thick red line) connecting 
them. It corresponds to a weakly coupled theory in the IR limit. The second UV complete 
trajectory runs towards large couplings, leading to a strongly coupled low energy theory 
with chiral symmetry breaking or conformality in the IR. The thin (magenta) trajectories 
emanating out of the UV stable hxed point correspond to the two irrelevant eigendirections 
of the hxed point. The other (thin magenta) trajectories illustrate that the Gaussian hxed 
point is an infrared hxed point for all couplings. Unlike in asymptotically free theories, it 
does not qualify as a UV hxed point. Finally, the thin (gray) lines show same trajectories 
which do not originate from a perturbative UV hxed point, but which nevertheless lead to 
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Figure 1: The phase diagram in the gauge-Yukawa-scalar subsector of couplings {ag,ay,ah) at 
NNLO accuracy with e = 0.05. Also shown are the UV and IR fixed points (dots), the UV safe 
trajectories (thick red line), a few distinguished trajectories (thin magenta lines) and sets of generic 
trajectories (thin gray lines). Arrows point towards the IR (see main text). 

a well-behaved weakly coupled theory at low energies. 

Figs. I^show the projection of Fig. [^onto the (a^, ay) and the (a^, ah) plane, respectively. 
Notice the near-perpendicular behaviour of trajectories close to the separatrix in the {ag, ay) 
plane (left panel). This is due to the fact that the velocity fidg of the RG flow of couplings 
along the separatrix is slower by an order in e <C 1. In the range of parameters considered 
here, the RG flow has a velocity of order e set by the Yukawa and scalar quartic couplings, 
yet the flow of the gauge coupling is smaller by a power in e. However, as soon as couplings 
are very close to the UV-IR connecting trajectory, the flow velocity becomes of order 
for all couplings, controlled by the gauge coupling. This effect is even more pronounced in 
the {ag,ah) plane (right panel), because the scalar and gauge /1-function are independent 
of each other at NNLO order. Gonsequently, the larger RG velocity of the scalar coupling 
becomes visible in the entire phase diagram, except along the separatrix. This effect plays 
also a role in the scalar subsector (ah, a^) for which two snapshots are shown in Fig. 
Here, the gauge and Yukawa couplings have been given values on the UV-IR connecting 
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Figure 2: The phase diagram of the gauge-Yukawa theory in the vicinity of the UV fixed point 
at NNLO accuracy with e = 0.05, projected onto the {ag,ay) plane (left panel) and the {ag,ah) 
plane (right panel). Shown are the asymptotically safe fixed point (UV), the Gaussian fixed point 
(G), and various sample trajectories. RG trajectories point towards the IR. 


separatrix. Specifically, we have taken Ug ~ 0.999 a* (left panel) and Ug ~ 0.397 a* (right 
panel), respectively, and ay determined via (2.15). Once again, in these projections, the RG 
flows in the scalar subsector have velocities of order e, whereas the flow velocity along the 
separatrix is of order and hence parametrically slower. Therefore the separatrix appears 
as a “pseudo” hxed point in this projection, indicated by a red dot. All other trajectories are 
attracted towards the separatrix. We also observe a second “pseudo” hxed point, indicated 
by a black dot. It relates to the trajectory which connects FP 2 with the Gaussian hxed 
point. However, along this trajectory the scalar held potential is unbounded from below 
and the theory is not considered to be physical. 


D. UV completion 


Trajectories which emanate out of the UV stable hxed point correspond to hnite theories 
at high energies. In the present case, the set of UV stable trajectories, sometimes denoted 
as the “UV critical surface”, is found to be one-dimensional [1]. The reason for this is a 
strictly dynamical one: the huctuations at the UV hxed point are such that only one linear 


combination of the four classically marginal couplings (2.2) becomes quantum-mechanically 


relevant. Similar dynamical constraints can also arise for theories with “complete asymptotic 
freedom” [8]. More specihcally, here, the one-dimensional critical surface leads to two UV 
safe trajectories. As indicated earlier, in the low-energy limit, they correspond to a weakly 
interacting theory of gluons, fermions and scalars with Gaussian scaling, and a strongly 
interacting theory with conhnement, chiral symmetry breaking, or conformality. 
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Figure 3: Projection of the phase diagram of the gauge-Yukawa theory onto the subspace of scalar 
couplings (NNLO with e = 0.05) with (ag,ay) taking values on the UV-IR connecting separatrix 
(left panel: Ug ^ 0.999 a* right panel; Og « 0.397ap. Red (black) dots indicate the trajectory 
which connect the physical fixed point FPi (unphysical fixed point FP 2 ) with the Gaussian fixed 
point. In the scalar subsystem where the RG flow is parametrically faster by 1/e, the separatrices 
appear as pseudo fixed points (see main text). 


The separatrix connects the UV hxed point with the Gaussian one. It coincides with the 
UV critical surface close to the hxed point P, and is characterised by global aspects of the 
phase diagram, namely the location of its hxed points and the eigendirections. Therefore, in 
general, identifying a separatrix requires global information about the phase diagram. For 
generic values of the parameter e, the separatrix can always be found numerically. For some 
considerations, it is useful to have explicit analytical expressions for the UV-IR connecting 
trajectory. Here, we explain how an analytical expression for the separatrix can be obtained 
locally to leading order in e 1 using the following observations. For small e 1, the four 
eigenvalues at the UV hxed point are given by 


^ 2 = %e + 0{e^) 

^ 3 = ^e + 0(e^) 

^4= ^V20 + 6V^e + O(e^). 


(2.13) 


We note that the irrelevant (positive) eigenvalues are of order e whereas the relevant (neg¬ 
ative) eigenvalue is of order e^. In consequence, the velocity of the RG evolution of 
couplings along the separatrix is oc e^, whereas the velocity of the RG evolution towards the 
separatrix is oc e. In the regime (2.8) we conclude that the approach towards the separatrix 
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is substantially faster than the evolution along the separatrix. For hxed this means 
that the corresponding value for the Yukawa coupling ay on the separatrix can be found by 
solving (3y = 0, to leading order in e. Following the same reasoning, for given (a^, ay) on the 
separatrix, the corresponding value for ah is then found from solving I3h = 0. Finally, a^ 
follows from solving I3y = 0 after inserting {ay, ah) on the separatrix. Following this logic, 
we are led to solve 

f3y = 0, 

f3h = 0, (2.14) 

/3v = 0 

for {ay, ah,ay) in terms of ag and to leading order in e. This leads to the following relations 
amongst couplings along the UV stable trajectories 


ah = ^{V^-l)ag , (2.15) 

ay = ^(^2oT6723 - 2^) ag . 


These relations for the separatrices are accurate in the limit e —)■ 0 (2.8). 
specify the RG running of the gauge coupling. 


It remains to 


E. Effective gauge coupling 


We now exploit the condition (2.14) to hnd the effective RG running of the gauge coupling 


along the separatrices. To simplify some of the subsequent expressions, we introduce a = ag 
throughout this section. Its RG flow is given by 


dta = -Ba^ + Ca^ + 0{a^). 


( 2 . 16 ) 


Below we will derive exact results valid for theories with (2.16) and generic B,C ^ 0. For 


the theory at hand, we concentrate on parameters in the regime B < 0 and G < 0 which 
take the specific values 


2 57 - 46e - 8e^ 
13+ e 


(2.17) 


They arise from (2.4) by expressing all a* {i = y,h,v) in terms of a, which are found by 


solving each of (2.14) to one-loop accuracy. Evidently, the effective gauge /3-function along 


the separatrix (|2.16|) displays three fixed points, a doubly-degenerated one at a* = 0, and 

a^, = B/C > 0 . 


an interacting one at 


( 2 . 18 ) 


Notice that in the absence of scalar helds and Yukawa interactions the coefficient C would 
read G = 25 > 0, and consequently the gauge coupling would not aquire an interacting UV 
hxed point in the physical domain a > 0. Including Yukawa interactions, the non-trivial 
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fixed point is perturbative as long as 0 < a* 1. In terms of (2.18), the solution to (2.16) 
is 

/ \ -B-at / \ / 

/i \ / a* — tto a \ /a* ct 

.ho 


a* 


a 


exp 


(2.19) 

«o y v^o 

Here, -B-a^ > 0, and a = a{fi) denotes the running gauge coupling along the separatrix, po 


the initial scale, and ao = a(/i = po) the gauge coupling at that scale. The solution (2.19) 
provides a closed expression for /r = One easily conhrms that (2.19) interpolates 

between the UV and IR hxed points. For a —)■ 0 the exponential term dominates, thereby 
showing that it corresponds to an IR hxed point with /r —)■ 0. In return, a —)■ a* corresponds 


to the UV limit where the algebraic term in (2.19) dominates, /r —)■ oo, and we recognize 
the appearance of the scaling exponent di = B ■ a^. Note that if ao > a* the Gaussian hxed 
point cannot be reached. Rather, RG trajectories then evolve towards strong coupling for 
all scales /i < /xq- In order to invert the solution p = fi{a) into a = a{fi) we introduce 


ct* 

IU = — 

a 


z = 


ho 

h 


-B-at 


a* \ /a* 

-1 exp-1 

ao J Vao 


( 2 . 20 ) 

( 2 . 21 ) 


For h/ho ranging between [0, cxd] and 0 < ao < we have that both ^ and W interpolate 
monotonously between [cx),0]. For 0 < a* < ao, we have that 2 ; becomes negative. Using 


these expressions, (2.19) turns into 


z = W exp W 


( 2 . 22 ) 


whose solution W = W[z] is given by Lambert’s W function. Uniqueness is observed 
provided that 2 ; > 0, which always holds true for the separatrix towards the Gaussian. In 
this regime the Lambert W function is monotonically increasing, polynomially so for small 
z and approaching 1112 ; from below for large z. For negative z, the Lambert W function is 
multivalued. The relevant branch is then the one with W[z] —?■ 0 for z —)• 0. 


In terms of (2.20) and (2.21), the fully resolved solution is hnally written as 


a(/i) = 


a* 


1 + IU(h) 


(2.23) 


where we introduced the shorthand notation IU(/i) = W[z{iJt)] with IU[ 2 ;] the Lambert 


function, and z{jj.) given in (2.21). Together with (2.15), the solution (2.23) fully describes 
the RG evolution of couplings along the separatrix to leading order in e. 


F. Characteristic energy scale and dimensional transmutation 

The dynamics of the theory towards low energies is characterised by an energy scale 
Ac where the RG flow displays a cross-over from UV hxed point scaling to either Gaussian 
scaling, or strong coupling, in the IR. The scale Ac arises dynamically through the mechanism 
of “dimensional transmutation”, and its role is similar to that of Aqcd for QGD. To identify 
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Figure 4: Renormalisation group running of all couplings along the UV-IR connecting separatrix 
(2.23), ( 2.15| ) in the NLO' approximation, with z{^) defined in ( |2.21 ). All couplings are normalised 
in units of a* (e = 0.05). The UV-IR crossover takes place at z{fj,) ^ 1. Although stays negative 
throughout, we notice that the sum of quartic scalar coupings ah + a^ stays positive, approching 
zero from above in the deep IR. 


the characteristic energy scale A^ we use(2.16) or (2.19) to introduce the effective scaling 
exponent 'd(a) = da(3g\^^^ which reads 


'd(a) = —2B a ( 1 — - — 

A Qijk 


(2.24) 


Close to the UV and IR hxed points, 'd reduces to the conventional scaling exponents 


^9(a^0) =-2-R-a^0+ 
'd{oL —)■ a*) = R ■ a* < 0 


(2.25) 


showing that the gauge coupling is marginally attractive at the Gaussian IR hxed point 
•Q —0+, and a relevant coupling at the interacting UV hxed point -& = B ■ oc —e^. In 


between the hxed points, (2.24) changes sign at 


a = - a*, 


(2.26) 


irrespective of the magnitude of B. The factor | arises because (2.24) vanishes at the point 
where the one-loop and the two-loop contribution, proportional to and a^, respectively. 


cancel out. In fact, all four /^-functions reach a maximum at (2.26) along the separatrix. At 


this point, and with decreasing RG momentum scale, or decreasing a, the RG how leaves the 
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basin of attraction of the UV fixed point and enters the domain of attraction of the IR hxed 
point. We define the characteristic energy scale as the scale where the UV-IR crossover 


happens. Using (2.19) and (2.26), can be expressed in terms of the renormalisation group 


scale parameter fi and a small deviation of a from the UV hxed point 5a = a* — a -C a* as 


Ac = c- iJ, 


1 - — 
a* 


(2.27) 


where we have suppressed terms of the order |5a|/a* 1. The proportionality constant is 

entirely hxed by the universal physics in the UV and, here, given by Ccross = exp[(ln2 — f)z/] 
with u = —l/'d(a*). Notice that the scale Ac is invariant under the renormalisation group. 


d , 

/U — Ac = 0 . 
dfi, 


(2.28) 


It arises from a dimensionless quantity, the small deviation |5a(/i)|/a* <C 1 from the UV 
stable hxed point at asymptotically large scales /r S> Ac through the phenomenon known as 
“dimensional transmutation”. For 6a < 0, and below the scale Ac, the theory is characterised 
by the infra-red free dynamics of weakly-interacting gluons, fermions and scalars, controlled 
by the Gaussian hxed point. A similar characteristic energy is linked to the UV hnite 
trajectory emanating from the UV hxed point towards strong coupling when /i —)■ 0 (5a > 0). 


The crossover away from the UV hxed point is characterised by the very same scale (2.27). 


Provided there is no strongly-coupled infrared hxed point we may estimate the scale where 


strong coupling sets in from ag(/i = Ac) ~ 1, hnding (2.27) with Cstrong = expi/. 


In the regime (2.8) , the results (2.27) for the characteristic energy scale and the crossover 


relation (2.26) are exact non-perturbatively, and to all loop orders.^ Once e is no longer 


taken to be asymptotically small, higher loop terms or even non-perturbative corrections 


will modify the relation (2.26), and, consequently, the numerical coefficient c in (2.27). 


For illustration of our results in this section, we display in Fig. the renormalisation 
group running of all couplings along the UV-IR connecting separatrix. In Fig. we have 


expressed all couplings as functions of 1/z, where z{jj.), dehned in (2.21) is related to the 


RG momentum scale. Using the expression (2.23) for the running gauge coupling, and the 
expressions (2.15) for the running Yukawa and scalar couplings to leading order in e, we 


observe that all coupling display a characteristic crossover from the interacting UV hxed 
point to the free IR hxed point. The crossover takes place at the RG invariant scale (2.27) 
where z ^ 1. We notice that a^, ay and ah stay positive throughout, whereas a^ stays 
negative. We also observe that the sum of the quartic scalar coupings ah + a^ stays positive 
for all scales. This result plays a central role for the stability of the theory to which we turn 
next. 


^ The result for the exact crossover relation (2.26) and the corresponding crossover scale Ac has already 
been used in m to explore asymptotically safe models of dark matter. 
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III. VACUUM STABILITY 

In the analysis of the renormalisation gronp flow it was tacitly assnmed that the vacunm 
of the scalar potential stays at the origin, such that all global symmetries are preserved 
along the flow. This hypothesis needs to be carefully scrutinised. In fact, we have already 
observed that the classical potential for one of the hxed points is unbounded from below 
and therefore the associated hxed point is not physical. We need to go beyond the classical 
stability for the other hxed point given that quantum corrections to the scalar potential can 
shift the minimum away from the origin. In perturbation theory this phenomenon is known 
as the Coleman-Weinberg mechanism m- This is therefore a crucial consistency check for 
the quantum stability of the vacuum at the UV stable hxed point and the associated critical 
how away from this hxed point. 


A. Classical moduli space 


We start by recalling the stability analysis for potentials of the form 

V = uTT(H^Hf+ v(TrB^ff)\ 


(3.1) 


where u and v are real couplings, and the x Np matrices H are invariant under U{Nf)lX 
U{NF)ji rotations. The main diherence to the original analysis of Paterson [19] is that the 
theory develops an interacting hxed point in the UV, where u and v should take their hxed 
point values. At this level, stability of the vacuum requires that we can, at most, have a set 
of flat directions along which we can hnd an inhnite number of vacua degenerate with the one 
at the origin. The set of hat directions - the classical moduli space - is parametrized by the 
matrices in held space with TrMMl = 1, V{M) = 0, and V (M -|- 5M) > 0 for all 5Mij. 
Since any classical held conhguration He can be diagonalized by a chiral U{Nf)l x U{Nf)r 


rotation, we can write He = diag(hi,..., h^p) without loss of generality. Consequently, (3.1) 
becomes 


Np / \ ' 

V = u^ht + vi^hU 

i=l \i=l / 


(3.2) 


If V vanishes for some He it will vanish for any multiple of He- It is thus sufficient to 
consider V on the hypersphere = 1- The stationary points of V on this hypersphere 

can be found by introducing a Lagrange multiplier —2A(J2i h‘f — 1) in (3.2) and solving the 
stationary equation: 


dV 

dhj 


Ahj \uh?j -|- n — a] =0 (no sum) . 


(3.3) 


The solutions are either hj = 0 or showing that all nonzero hi must be equal at 

an extremum. From the hypersphere constraint it follows that for n nonzero elements hi, 
each of these must be related to n as flf = 1/n. Evidently n may take values between 1 
and Nf- Consequently, at an extremum the potential reads Wxtr = v + u/n. The value of 
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fixed point 

single trace coupling 

stability condition 

stability fullfilled 

FPi 

< >0 

al + al>0 

yes 

FP2 

<>o 

al + al>0 

no 

FP3 

< <0 

al -b a.v{fi)/Np > 0 

infinite Np : no 

finite Np : possibly 


Table II: Classical moduli space. Full stability of the scalar potential at infinite Np is achieved 
at the UV fixed point (2.9) of the gauge-Yukawa theory. Notice that a second domain of stability, 
related to the partial fixed point FP 3 , opens up for finite Np and sufficiently large a^,. 


n is determined by requiring the extremum to be a minimum. This can be realized in two 
different manners: Provided that m > 0, then Vextr is minimal for n = Np, and consequently 
Mij = 6ij/\/Np. Alternatively, provided that u < 0 , then Vextr is minimal for n = 1 , and 
consequently = diag(l,0,... ,0). We conclude that for V to be bounded from below, 
the couplings are constraint to the parameter region {u > 0 A n - 1 - u/Np > 0 for oc 6ij}, 
or to the region {M<0Ar’-l-M>0 for Mij oc 5^1 }• Flat directions are obtained if, in either 
case, the second inequality becomes saturated. 


B. Vacuum stability at UV fixed points 


Using the hndings of Sec. |III A[ we now turn to the UV hxed points detected in our 


model. In the conventions introduced in (2.2), and depending on the sign of the single-trace 


scalar coupling, the viable domains correspond to either of 

ah > 0 and an + a^ > Q , (3.4) 

or ah < 0 and ah + a^/Np > 0. (3.5) 

The stability analysis for all hxed points is summarised in Tab. |TTj Both FPi and FP 2 have 
ah > 0, but the condition (3.4) only holds true for FPi. At the partial hxed point FP 3 , we 


observe ah < 0, but > 0 displays a Landau pole towards high energies. However, in the 


strict inhnite Np limit, a^ drops out of (3.5) and consequently, (3.5) cannot be satished for 
any ah- We conclude that of all hxed points within the perturbative regime, it is FPi (2.9) 


which corresponds to a stable scalar held potential, whose global minimum is located at the 
origin in held space (see Tab. 0 , Furthermore, the potential has no hat directions provided 
that e > 0. At hnite Np, however, additional stable solution of the type (3.5) may arise 
in the regime of large a-u- Presently this cannot be decided based on the approximations 
adopted here. 


C. Quantum moduli space 

We now investigate the quantum moduli space through the explicit computation of the 
quantum ehective potential a la Coleman-Weinberg [niED]. We are particularly interested 
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in the UV safe trajectory which connects the UV fixed point with the Gaussian IR fixed 
point. From the moduli-space analysis in the previous subsection we can assume, without 
loss of generality, that Hij —)■ (pcSij where the constant classical field 0 c takes the role of 
the expectation value of the renormalised quantum field. The quantum effective potential 
thereby becomes a function of the renormlasied field, the renormalisation group scale, and 
the running couplings. As such, it obeys an exact renormalisation group equation 


d d d 

-7(0^) (pc— + ;i0,Oj) = 0. 


(3.6) 


duo d(j)c ^ doti 

Similar exact renormalisation group equations can be written down non-perturbatively for 
general effective potentials [211122]. Here, the functions f3i are given in (2.4) - (2.7), and the 
function 7 = —^dlnZ/dln/i denotes the scalar field anomalous dimension given by 

7 = CTy — § (^ + e) Q!y + ^ayag + 2a\ + 0{a^) (3.7) 

up to the two-loop order in perturbation theory [T]. To solve the exact renormalisation 
group equation (3.6), we exploit that the potential is a pure quartic at the fixed point. By 
dimensional analysis the effective potential can then be written as 


heff (0CJ hO) ^es(^4^c/ h'O) ^i) ' 0c 


(3.8) 


in the absence of further mass scales. The dimensionless prefactor Aeff(0c//io, can only 
be a function of dimensionless parameters such as the dimensionless couplings, and the 
background field in units of the renormalization scale /iq. It obeys the exact renormalisation 
group equation 



+ 47(«j) - 

i 



'^eff(0c) 0 1 


(3.9) 


whereby we have traded the pcl^-dependence for the 0 cl 0 -dependence. We have also intro¬ 
duced the new functions 




7(tti) 


Pi 

1 -h 7 (a 0 
7(«i) 

1 -h 7(a0 ’ 


(3.10) 


which allow us to translate information about the renormalization group 0 -functions of 
couplings into the field-dependence of the potential. Integrating the renormalization group 
equation for Aefr( 0 c) with a suitable boundary condition Aefr( 0 c = /^o) is then equivalent to 
a resummation of logarithms in the field. Technically, this is done by defining new running 
couplings ai through the differential equations 


^ ai{4>c ^ Mo) ^ a* 


(3.11) 
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Figure 5: Stability of the quantum vacuum and absence of spontaneous symmetry breaking: shown 
is the Coleman-Weinberg effective potential 14fr(<^c) along the UV finite trajectory connecting the 
UV stable fixed point and the Gaussian IR fixed point, as a function of the constant classical field 
0c (e = 0.05). The potential is normalised to the value of the classical fixed point potential at 
0c ~ /^o- 


This allows us to write a general solution of (3.9) in terms of the couplings a* which solve 
as 


Aeff((/>c) = A(0c) exp -4 


' Mo 


du . 
— 7(r) 

/i 


(3.12) 


where A (0c) = A [a* (0c)] can be an arbitrary functional of the functions whose scale- 
evolution is determined via (|3Tg. Its field-dependence solely arises implicitly through 
the dependence of di on ln(0c/po)- The functional form for A(0c) can be hxed through a 
matching with the hxed order perturbative expansion. To one-loop accuracy, and recalling 
the normalisation convention (2.2), the function A(0c) reduces to the sum of the trace and 
double-trace scalar couplings, 




(3.13) 


where the substitution ai{(j)c) -d di[(l)c) ensures that the leading logarithms are resummed. 
The validity of conventional perturbation theory requires the smallness of both \ai\ 1 
and Ittj ln(0c//io)| ^ 1. The resummation reduces the two constraints to a single one, 
|tti(0c)| ^ 1- Consequently, the regime of small helds can now be probed even though the 
logarithms ln(0/po) becomes large in the limit of small helds.^ 


^ An example showing that the fixed order effective potential without resummation is unable to capture 
the physical properties of the potential near fixed points is given in [23] . 
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D. Analytical Coleman-Weinberg potential 


We now derive analytical expressions for the resummed quantum effective potential to 
NLO' order in the approximation, the leading non-trivial order in e where the RG flow of 
scalar couplings becomes relevant. We can substitute 7 by 7 and suppress the subleading 7 - 


corrections in (3.10). In this limit, the explicit solution (2.23) for the running gauge coupling 


is valid, and the running of the scalar couplings along the separatrix strictly follows the 


running of the gauge coupling via (2.14). To hnd an explicit expression for the exponential 


term in (3.12), we recall from (3.7) that 7 = at one-loop. By the virtue of (2.14), 


this leads to j = D ag with H ^ on the trajectory of interest ( 2.14[ ). The remaining 
integration is then performed analytically using the explicit form for the running gauge 


coupling (2.23), together with (2.20) and (2.21). The classical and the resummed quantum 


effective potentials then take the explicit form 


Ki(0c) = A, 0^ 


-AD/B 


(3.14) 


where the positive quartic coupling is given by A* = e ^^(\/20 + 6 a/^ — — 1 ) at 


the hxed point, see (2.9), and the exponent —AD/B = 18/(13 ■ e) > 0 is positive and 


parametrically large in the regime of interest (2.8). For the computation of We we have 


adopted the UV boundary condition where Aeff(</>c = fJ-o) is very close to A* to ensure that 
couplings are in the close vicinity of the UV hxed point. This implies that W{fiQ) -C 1. 
The effective potential then reduces to the classical potential at the hxed point in the limit 
(/ —)■ /To- Corrections due to the resummation of logarithms arise once (/c < ho, and they can 
become sizeable once 0c "C ho- The tendency of these corrections can be understood from 
(3.14) as follows. With decreasing 0c < ho, we have that lU(0c)/kF(ho) ~ (ho/0c)~'^“*- 
This leads to a tiny reduction cx (1 — lU(ho)(ho/0c)~^“*) of Ves over 101 due to the running 
of the gauge coupling, and to an enhancement oc (ho/ 0 c)^^®'^^^*'^^^ due to the anomalous 
dimension of the scalar held. Evidently, the enhancement wins, meaning that the resummed 
logarithmic corrections increase Ves over 10i with decreasing 0c/ho- We also note that the 
shape of the potential is entirely dominated by the vicinity to the UV hxed point. It is only 
for held values 0c of the order of the scale Ac (2.27) and below where the Gaussian hxed 


point takes over the control of the resummation of logarithms. 

For want of completeness, we now investigate the corrections to the ehective potential 


(3.14) due to the additional resummation of the scalar held anomalous dimension in (3.10), 


(3.11). Integrating (3.11) for using the two-loop result for and the one-loop expression 


for 7 leads to the explicit solution 


a* 


a(0c) = 


1 + lU(0c) 


(3.15) 


for the gauge coupling, with kF(0c) = lU[2;(0c)], where W is the the Lambert function and 
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z{fi) given by the expression 


z = 


ho 

h 


—B-at/{l+D-a,) 


a* 

— _ 1 
ttO 


exp 


^ _ 1 
bo 


(3.16) 


Notice that in (2.21) is identical to 2 : in (3.16) after the replacement Ba^ —)■ Ba^/{1+D a*). 


The effect of the additional resnmmation, therefore, is that the rnnning of the gange conpling 
dg{4>c) from the UV to the IR is mildly accelerated over the rnnning of ag(0c) at the same 
scale 0c. Since both conplings are monotonically decreasing fnnctions along the UV-IR 
connecting trajectory, this implies that bg(0c) < «g(0c)- Eqnality holds only at the IR and 
UV hxed points, corresponding to the limits 0c 0 and 0c —)■ ho respectively. Fnrthermore, 


the relations (2.15) continne to hold true for the couplings a* because all functions 0, follow 


from 0j through a rescaling by one and the same factor. For this reason, A (0c) takes the 
same functional form given in (3.13), with ah and Uy expressed as functions of dig in (3.15), 


(3.16). Similarly, the exponential factor in (3.12) can also be integrated analytically for 7 


given in (3.10). The effective potential then takes the form 

Ucl(0c) 


Ueff(0c) = — ^ = 

l + lU(0c) VlV(ho) 


lU(0c)^ 


1 + Z) a* + W {(j)c) 
1 + Z) a* + W (ho) 


-4 


B l + D-n* 


(3.17) 


Comparing (3.17) to (3.14), and also (2.21) with (3.16), we observe that the additional 


resnmmation has induced corrections of the order of D a* 
corrections are parametrically suppressed by e in the regime (| 2 . 8 |), and, additionally, by a 


Since D a^ = these 


factor of Quantitatively, (3.17) and (3.14) are essentially indistinguishable for helds 0c 


larger than the characteristic scale Ac (2.27) where the underlying RG flow crosses over from 
UV to IR scaling. Once 0c is of the order of Ac and below, meaning W,W ^ h and above. 


the potential (3.17) is enhanced over (3.14) starting in the 5-10% regime (for e = 0.05). We 
conclude that the additional corrections due to 7 —)■ 7 and dt in (3.10), (3.11) are 

minute and quantitatively negligible unless 0c Ac, and we can stick to the result (3.14) 
for all technical purposes. 


E. Coleman-Weinberg potential and strong coupling 


As indicated earlier, two UV complete trajectories emanate from the UV stable hxed 
point. The hrst UV complete trajectory, towards weak coupling, connects the UV hxed 
point with the Gaussian IR hxed point. The cross-over to Gaussian behaviour takes place 
at the RG invariant scale (2.27). The second UV complete trajectory runs towards a regime 


of strong coupling in the IR. Again, the cross-over towards strong coupling takes places 
at a scale Ac as in (2.27). Interestingly, the results of the previous subsection for the 
Goleman-Weinberg potential are applicable for either of these trajectories. The main dif¬ 
ference between them is that the gauge coupling decreases in the former case, but grows in 


the latter. As we have discussed in Sec, HID close to the UV hxed point the running of the 
gauge coupling is too slow to leave an impact on the Goleman-Weinberg potential. Rather, 
it is the sign of the scalar held anomalous dimension which controls the logarithmic quan¬ 
tum corrections. To leading order, and for helds 0c larger than the characteristic scale Ac, 
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Figure 6: Enhancement of the quantum effective potential over the classical one: shown are the 
resummed quantum effective potentials 14fr(<^c) at different orders in the approximation (NLOh 
dashed line, NNLO: full line) for e = 0.05, and normalised to the classical potential lci((^c) at 
the UV fixed point, see (3.14). At leading logarithmic accuracy, the enhancement 14fj^NNLo('?^c) > 
Kff,NLO'('/’c) > Vci{(j)c) is due to ^ 7nlo ^ 0 (see main text). 


these corrections are largely independent of whether the gange conpling grows or decreases. 
Hence, with initial conditions ao at /iq on either side of, bnt very close to, a*, the UV hxed 
point controls the logarithmic corrections and the corresponding effective potentials become 
indistingnishable from each other. On the other hand, for much smaller helds, say of the 
order of 0 c ~ Ac, subleading differences between the potentials for weak and strong gauge 
coupling start to become visible. In this regime, perturbation theory is no longer a good 
approximation along the strong-coupling trajectory. For either of these trajectories, we con¬ 
clude that quantum effects neither destabilise the symmetric vacuum, nor do they lead to 
the spontaneous breaking of symmetry. 


F. Numerical Coleman-Weinberg potential 


We now turn to a full numerical evaluation of the quantum effective potential at NNLO 
accuracy, corresponding to (3,2,l)-loop accuracy in the gauge, Yukawa, and scalar sector. 


respectively, using the /3-functions as given in (2.4) 


(2.7). Unlike the preceeding approx¬ 
imation, the RG flows at these orders obey Weyl consistency conditions. Concretely, for 
e = 0.05, we determine the UV hxed point and the UV-IR connecting trajectory numeri¬ 


cally for all couplings. These are then exploited to hnd the quantum effective potential (3.8) 


by numerical integration of (3.11) and (3.12) 


Our hndings are shown in Fig. [^in comparison with the analytical results at NLO' order 
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(3.14) and (3.17). As can be read off from Fig. with decreasing 0c/hO) the effective 
potential Ve^ at NNLO (fnll line) is enhanced over the NLO' resnlt (long dashed line). The 
red short dashed line indicates the result without resummation of logarithms. The slight 
enhancement of the NNLO result over the NLO' result is understood as follows. The main 
effect relates to the scalar anomalous dimension 7 = whose value 7nnlo is enhanced over 
7nlo due to a slight shift in the hxed point value of ay from NLO' to NNLO. At e = 0.05, the 


difference is approximately 20%. In view of (3.12), (3.13), this leads to the relative change 


14 


eff.NNLO 


14 


efr,NLO' 


_L //, \ 4(7nnlo-7nlo) 

+ a^jNNLO / ho \ 

+ a*)NLO' \4>cJ 


(3.18) 


to leading logarithmic accuracy in (/)c/ho- The estimate (3.18) is in very good quantitative 


agreement with the hndings of Fig. A secondary effect relates to the running of the gauge 
coupling which at NNLO is slightly slower that at NLO'. This is evidenced in Fig. and by 
the fact that the relevant scaling exponent has changed by roughly 10% from hi ~ —0.0015 
at NLO' to hi ~ —0.0014 at NNLO. In principle this would lead to a small reduction of 


the ratio (3.18), which, however, is exponentially small and hence of no relevance for the 


potential in the region considered here. We can therefore conclude that the enhancement 


14l(0c) < Kff.NLO' (0c) < 14fr,NNLo(0c) (3.19) 

of the effective potential for (pc < Po to leading logarithmic accuracy is largely due to the 
scalar anomalous dimension whose hxed point value 0 < 7nlo < 7nnlo changes from order 
to order in the approximation. Evidently, in this theory, quantum effects neither destabilise 
the symmetric vacuum, nor do they lead to the spontaneous breaking of symmetry. 

Finally, we discuss our results in the light of perturbation theory and Weyl consistency 
conditions. The interacting UV hxed point arises for the hrst time at NLO’ accuracy, 
see Tab. |T} Stability of the vacuum is observed, persisting to hnite values of e. Stability 
persists beyond NLO’ accuracy, where non-universal parameters, such as the gauge three- 
loop coefficients or the Yukawa two-loop coefficients, come into play. This strengthens the 
view that the hxed point exists for hnitely many helds. Moreover, we observe that the 
stability of the vacuum is compatible with Weyl consistency conditions. The latter arise 
from formal considerations of path integrals on curved backgrounds, and relate scheme- 
dependent orders in perturbation theory, e.g. the NLO and the NNLO approximation. Weyl 
consistency conditions are prima facie insensitive to physical observables such as minima of 
ehective potentials. In this light it is noteworthy that Weyl consistency is compatible with 
vacuum stability. 


IV. CONCLUSIONS 

Recently, novel classes of fundamental four-dimensional quantum held theories with non- 
Abelian gauge helds, fermions, and scalars have been discovered whose high-energy be¬ 
haviour is asymptotically safe, controlled by an exact interacting UV hxed point [T]. Renor- 
malisation group trajectories emanating from the hxed point relate to well-dehned, hnite. 
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Figure 7: Shown is the gauge /3-function along the UV-safe trajectory at NLO (dashed line) and 
NNLO (full line) accuracy, with e = 0.05. The main NNLO effect consists in a slight reduction of 
the /3-function once couplings are away from their fixed points. 


and predictive local quantum field theories at all energies, despite of the fact that asymptotic 
freedom is absent. The hxed point occurs parametrically close to the Gaussian and admits 
rigorous control within perturbation theory. 

We have extended the study of pQ to establish that the vacuum of UV safe gauge-Yukawa 
theories is stable, classically and quantum-mechanically, even though asymptotic freedom 
is absent. We also found that the main quantum corrections to the effective potential arise 
due to the anomalous dimension of the scalars. Unlike in asymptotically free theories, here, 
the scalar anomalous dimension takes a non-vanishing value even at highest energies. The 
renormalisation group running of couplings away from the hxed point is a subleading effect 
for the effective potential, provided held values remain large compared to the characteristic 
energy scale Ac (2.27) of the theory. The absence of classically hat directions of the hxed 
point potential thus entails quantum stability. Owing to the perturbative nature of the 
hxed point, we also determined the crossover relation (2.26) exactly, including analytical 
expressions for all running couplings along the UV safe trajectories. 

We have limited our investigation to massless theories in the Veneziano limit where the 
number of helds is very large and UV interactions are weak. Continuity in the number 
of helds indicates that the vacuum remains stable even for hnitely many helds as long 
as perturbation theory remains a good approximation. At strong coupling, functional 
renormalisation can be used to access the Coleman-Weinberg potential non-perturbatively 
[21]. Our proof of vacuum stability can straightforwardly be exported to other gauge 
theories with interacting UV hxed points. 
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